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Abstract. The RR Lyrae instability strip (IS) in NGC 1851 is investigated, and a model is 
presented which reproduces the pulsational properties of the RR Lyrae population. In our 
model, a stellar component within the IS that displays minor helium variations (y~0.248- 
0.270) is able to reproduce the observed periods and amplitudes of the RR Lyrae variables, 
as well as the frequency of fundamental and first-overtone RR Lyrae variables. The RR 
Lyrae variables therefore may belong to an He-enriched second generation of stars. The 
RR Lyrae variables with a slightly enhanced helium (y~0.270-0.280) have longer periods 
at a given amplitude, as is seen with Oosterhoff II (OoII) RR Lyrae variables, whereas 
ly*-} , the RR Lyrae variables with y~0.248-0.270 have shorter periods, exhibiting properties of 

■ Oosterhoff I (Ool) variables. 
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1. Introduction have different C+N+O contents (e.g., Cassisi 

et al. 2008; Ventura et al. 2009). The horizontal 
The globular cluster (GC) NGC 1851 hosts branch (HB) of NGC 1851 also hosts two dis . 

two distinct subgiant branches (SGBs) (Milone tinct popu i ations? name l y a prominent red HB 

et al. 2008), prompting much discussion as clump and a blue ta0 Finallyi differ e n t popu- 

to the formation history of this cluster. The lations have also been reve aled in the red giant 

X ! split between the bright SGB (SGBb) and faint branch (RGB) (e g ^ Gmndahl et aL i 999 - Han 

5— i . SGB (SGBf) can be explained by two sub- et a j 2009) 

populations that differ in age by about 1 Gyr 

(e.g., Milone et al. 2008; Gratton et al. 2012), xhe HB of g i obu i ar cms ters (GCs) is a 

or by two SGBs that are nearly coeval but va i uab i e stellar component which can be used 

to investigate the formation and evolution of 

Send offprint requests to: A. Kunder GCs (e.g. Gratton et al. 2010; Dotter et al. 
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2010). Indeed, there have been previous stud- 
ies to model the HB of the NGC 1851 to ob- 
tain scenarios of the formation this cluster's bi- 
modal horizontal branch (Salaris et al. 2008; 
Gratton et al. 2012). Our study concentrates on 
the portion of the HB inside the instability strip 
(IS), by discussing the case of the RR Lyrae 
properties in NGC 1851. 

2. RR Lyrae Observations 

Because of the compact nature NGC 1851, it 
is likely that the RR Lyrae star sample in the 
crowded core in incomplete (e.g. Sumerel et al. 
2004; Downes et al. 2004). Fortunately at dis- 
tances greater than 40" from the cluster cen- 
ter, individual stars can be relatively easily 
resolved, and there is no indication that the 
RR Lyrae sample is incomplete in this region. 
Therefore we limit our sample of RR Lyrae 
stars with which to compare our HB models 
to this outer region. 

The outer RR Lyrae variables were stud- 
ied by Walker (1998) using the BVI pass- 
bands, and robust mean magnitudes, periods 
and amplitudes are available. The periods and 
V-amplitudes of our sample of fundamental 
mode RR Lyrae (RR0) variables are shown 
in Figure 1, and the period-amplitude relation 
of typical Ool and OoII-type systems is over- 
plotted. We note that although many of the RR 
Lyrae stars in NGC 1851 have periods and am- 
plitudes that cause them to fall near the Ool PA 
relation, there are a number of stars following 
the OoII PA relation. 

In comparison, Figure 1 shows the peri- 
ods and V-amplitudes of 1097 RR0 Lyrae vari- 
ables in 40 Galactic globular clusters (Kunder 
et al. 2013). This sample of RR0 Lyrae 
variables is divided by their position in the 
period-amplitude plane following the lines 
that Clement & Shelton (1999) derived for 
Oosterhoff I and Oosterhoff II RR0 stars (see 
Figure 1). Those RR Lyrae variable falling 
closest to the Ool line are designated as Ool- 
type RR Lyrae stars, and those falling closest 
to the OoII line are the OoII-type variables. 

We define an Oosterhoff ratio for each GC, 
which is simply the number of Ool-type RR0 
Lyrae stars compared to the total number of 
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Fig. 1. Top: Period-amplitude diagram for our sam- 
ple of RR0 Lyrae variables in NGC 1851. Bottom: 
Period-amplitude relation for 1097 RR0 Lyrae vari- 
ables in 40 Galactic GCs. Our division of Ool and 
OoII-type RR Lyrae variables are shown by dark 
and light circles, respectively. The lines derived by 
Clement & Shelton (1999) for Ool and OoII RR0 
are overplotted. 

RR0 Lyrae stars in the GC, OoIrro/Toirro- 
Figure 2 shows the histogram of the Oosterhoff 
ratio of the GCs in our sample. Our defined 
Oosterhoff ratio splits the Milky Way GCs into 
two groups; the Ool-type clusters have RR 
Lyrae variables with shorter periods for a given 
amplitude and hence have larger Oosterhoff ra- 
tios (with respect to the OoII-type clusters). 
Further, there is an absence of clusters falling 
in the "gap". We therefore believe that our 
Oosterhoff ratio is useful to distinguish be- 
tween Ool- and OoII-type GCs. 

3. Synthetic HB Modeling 

Our synthetic HB calculations are carried out 
to provide a simple and attractive explana- 
tion for the cluster HB and IS morphology, 
keeping the number of free parameters to a 
minimum, yet still reproducing the RR Lyrae 
properties that make this cluster stand out as 
having an unusual Oosterhoff type. The HB 
evolutionary tracks used are from the BaSTI 
stellar library (Pietrinfemi et al. 2004, 2006, 
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Fig. 2. A histogram of the ratio of Ool stars in 
the Milky Way GCs. The Oosterhoff ratio of 
the NGC 1851 RRO Lyrae stars is highlighted. 



2009). The horizontal part of the HB, which 
includes the RR Lyrae instability strip, makes 
up ~10% of the cluster stellar sample brighter 
than V — 21 mag; this is the component our 
study focuses on modelling. To be considered 
a match to the observations, the synthetic HB 
model must agree with four specific observa- 
tional constraints. The first is that the num- 
ber ratio of the progeny of the SGBb to the 
progeny of the SGBf subpopulations has to sat- 
isfy the 70:30 ratio observed along the SGB 
(Milone et al. 2009). The second is that the 
(B :V :R) (blue :variable :red HB) ratio agrees 
with the observed (B :V :R) = (33 + 8:10 + 
5 :56 ±11) (in line with the results by Catelan 
et al. 1998; Saviane et al. 1998). The third re- 
quirement is that the theoretical periods and 
amplitudes agree with those observed, and the 
fourth that the frequency of first-overtone RR 
Lyrae stars (RR1) to RRO variables is approx- 
imately NJN tot ~ 0.35 (Walker 1998; Downes 
et al. 2004). 

The four HB components described by 
Gratton et al. (2012) are used as a starting 
points for our calculations. As in Gratton et al. 
(2012), The only parameter we vary is the ini- 
tial He abundance of the HB progenitors, keep- 
ing the same total RGB mass loss for all com- 



ponents. That the mass loss is the same for all 
components is not something we can prove, 
and we stress that this is our working assump- 
tion, and our results follow from that. 

Objects from our synthetic HB that fall 
within the observed IS from Walker (1998) 
are considered RR Lyrae variables. Gratton 
et al. (2012) model the horizontal part of the 
HB component with a constant He abundance 
7=0.265. We find a much better fit to the RR 
Lyrae properties, however, by assuming a con- 
tinuous He distribution between 7=0.248 and 
0.280 (see Figure 3). In our simulation, the 
mean He abundance in the IS is < 7 >=0.271, 
close to the constant abundance 7=0.265 em- 
ployed by Gratton et al. (2012) for this compo- 
nent, and the mean mass is < M > = 0.634 Mq. 

Figure 3 shows the periods of the RR Lyrae 
variables in our simulation as compared to the 
observed periods, where the theoretical periods 
are calculated for all HB stars falling within 
the observed IS using the Di Criscienzo et al. 
(2004) RR Lyrae pulsation models. The RR1 
variables are fundamentalized via log Po ~ 
logP+0.127, where Pq is the fundamental mode 
period. A Kolmogorov-Smirnov (KS) test be- 
tween the observed periods and the synthetic 
RR Lyrae periods returns a probability P=0.86, 
well above the default threshold P t h=0.05 be- 
low which one rejects the null hypothesis. On 
this basis, we find that the synthetic periods 
from our simulated HB and the observed pe- 
riods agree well with each other. 

The ratio of first overtone to total RR Lyrae 
variables in our simulation is N\ /N tot ~0.30, in 
agreement with that observed. In contrast, sim- 
ulations using a constant helium for the portion 
of the HB containing the IS (as in Gratton et al. 
2012) do not fit the constraints given by the 
NGC 1851 RR Lyrae variables as well. For ex- 
ample, adopting 7=0.265 results in an Ni/N tot 
= 0.1 1 (versus the observed Ni /N tot = 0.28). 

4. Conclusions 

Modern theoretical synthetic HB models are 
used to investigate the RR Lyrae properties of 
NGC 1851. Both the RR Lyrae period distri- 
bution and the number ratio of first overtone 
RR Lyrae to total RR Lyrae stars, Ni /N tot , pro- 
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Fig. 3. Top: A comparison between the ob- 
served NGC 1851 RR Lyrae periods versus the 
theoretical periods with a constant Y=0.265 
(Gratton et al. 2012). Bottom: The results when 
using a He abundance with minor helium varia- 
tions (Y~0. 248-0. 280) versus the observed pe- 
riods. 

vide constraints pertaining to the component of 
the HB containing the IS. It is straightforward 
to reproduce the observed distribution of RR 
Lyrae stars inside the instability strip with mi- 
nor He variations (7~0.248-0.280). 

In general, the RR Lyrae variables with 
Y < 0.27 fall in the Ool area of the PA dia- 
gram, whereas the RR Lyrae variables with Y 
> 0.27 fall close to the OoII line. Assuming 
that the period-amplitude diagram can be ef- 
fectively used to classify RR Lyrae stars into 
an Oosterhoff type, this means that He and 
Oosterhoff type are correlated in this clus- 
ter. This is not completely unexpected, as 
an increase in He makes RR Lyrae variables 
brighter and, as a consequence, higher he- 
lium abundance makes their pulsational period 
longer at fixed temperature/color (Bono et al. 
1997). 
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